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Abstract 

We show that non-Hcrmitian and nearest-neighbor-interacting perturbations to the Fritzsch 
textures of lepton and quark mass matrices can make both of them fit current experimental 
data very well. In particular, we obtain — 45° for the atmospheric neutrino mixing angle 
and predict ~ 3° to 6° for the smallest neutrino mixing angle when the perturbations in 
the lepton sector are at the 20% level. The same level of perturbations is required in the quark 
sector, where the Jarlskog invariant of CP violation is about 3.7 x 10^^. In comparison, the 
strength of leptonic CP violation is possible to reach about 1.5 x 10^^ in neutrino oscillations. 
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1 Introduction 



The flavor sector in the standard model (SM) of electroweak interactions has been puzzling because 
it involves most of the free parameters of the model itself. Although the values of six quark masses 
and those of three angles and one CP-violating phase of the 3x3 Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix V [1] are all known to a good degree of accuracy [2], it remains very 
difficult to understand why they have the observed mass spectrum and flavor mixing pattern. The 
situation in the lepton sector is even worse: not only the absolute mass scale of three neutrinos but 
also the smallest angle and CP-violating phases of the 3x3 Maki-Nakagawa-Sakata-Pontecorvo 
(MNSP) lepton mixing matrix U |3J are still unknown. To resolve the flavor problem in the SM 
one has to gain an insight into the flavor structures and possible flavor symmetries behind them. 
In spite of many attempts in this direction, a successful (unique and predictive) flavor theory has 
not been achieved. Most of the present-day studies on the flavor structures of quarks and leptons 
are more or less phenomenological and among them the texture- zero approach [5] has proved 
to be very useful to establish some simple and testable relations between the mass ratios of quarks 
or leptons and their corresponding flavor mixing angles. 

In the three-family framework the Pritzsch texture of quark mass matrices [6] 



Mf ) = A* i? J , (1) 




where a = u (up) or d (down), has attracted a lot of interest since it was proposed in 1978. It 
belongs to the more generic nearest-neighbor-interaction (NNI) form of quark mass matrices, 



Mr'^ = A'^ b\ . (2) 




The NNI form can always be obtained from an arbitrary form of My and via a proper choice of 
the flavor basis in the SM [7] . So the Pritzsch texture is actually a NNI texture with the additional 
assumption of the Hermiticity conditions A'^ = ^* and B'^ = 5*. In view of the problem that 
Mu and cannot simultaneously give rise to a sufficiently large value of the top-quark mass 

and a sufficiently small value of the CKM matrix element | V^^ | , one has to abandon either the 
NNI feature of quark mass matrices [8] or their Hermiticity [9], either partly or completely 0. It 
is always possible to numerically determine the departures of realistic and from and 
by using current experimental data. If such departures are not very significant, they can be 
regarded as small perturbations and then be treated in an analytical way so that their effects on 
the CKM matrix elements will become more transparent. 

It also makes sense to consider non-Hermitian and nearest-neighbor-interacting perturbations 

(V) (V) 

to the Fritzsch-type lepton mass matrices and M,) . Although the latter can fit current 

experimental data [10], it is very difficult to obtain — 45° or equivalently the maximal or 
nearly maximal atmospheric neutrino mixing. One possible way out is to introduce the seesaw 



^One may certainly abandon both the NNI and Hermiticity conditions by following a different starting point of 
view (e.g., the triangular form of and M^i) to investigate quark mass matrices and their consequences on flavor 
mixing and CP violation 4 . But we do not focus on this possibility in the present work. 
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mechanism to the neutrino sector [llj, in which the Dirac neutrino mass matrix Mj-, takes the 
Pritzsch texture as does but the heavy Majorana neutrino mass matrix AI^ is (approximately) 
proportional to the identity matrix. Here we follow a different way at the electroweak scale. We 

(V) (V) 

shall introduce non-Hermitian perturbations to both ' and Mj) ' so that the resultant charged- 
lepton and neutrino mass matrices can agree with current neutrino oscillation data to a better 
degree of accuracy [^. Such a parallel study of lepton and quark mass matrices of approximate 
Fritzsch textures is useful to reveal the similarities and differences between the lepton and quark 
sectors, and it should also be helpful for building a unified flavor model of leptons and quarks. 

Let us point out that the present work is different from a recent one done by Branco et al [12] in 
the following three aspects. (1) They have only considered quark mass matrices of the NNI form, 
whereas we are discussing both lepton and quark mass matrices of the same NNI form and giving 
a stronger emphasis to the lepton sector. In particular, we are concerned about an interpretation 
of the observed ^23 — ^5° ^^^^ a prediction of nonzero ^^^3 for the MNSP matrix based on the 
NNI texture of lepton mass matrices at the electroweak scale. (2) The analytical approximations 
made in our perturbative calculation are valid to a better degree of accuracy and thus allow one 
to see the difference between the contribution of fermion mass ratios and that of perturbation 
parameters (which signify a departure of the NNI texture from Hermiticity) to the flavor mixing 
matrix elements in a clearer way. (3) Our numerical analysis is more comprehensive than the one 
done in Ref. [12j, and it shows that current experimental data require the non-Hermitian effects 
to be at the 20% level in both lepton and quark sectors. This observation is expected to be useful 
for model building, especially when leptons and quarks are discussed in a unified flavor picture. 

The remaining part of this paper is organized as follows. In section 2 we do a perturbative 
calculation to reveal the salient features of non-Hermitian corrections to the Fritzsch textures of 
lepton and quark mass matrices. Section 3 is devoted to a numerical illustration of the constrained 
parameter space at a reasonable level of perturbations in the lepton and quark sectors. A brief 
summary, together with some further discussions, is given in section 4. 



2 A perturbative calculation 



Without loss of generality, the NNI mass matrix m'^^'^'^ (for a = u, d, / or v) in Eq. (2) can always 



be decomposed into M, 
matrices, and 



(NNI) 



P^m2^^^^ P'^, where and P'^ are two independent diagonal phase 



M(NNI) 



/O 




b' 



(3) 



is real. After the bi-unitary transformation OaMa^^^^O^ = M„ = Diag{A", , A3 } with Af (for 



i = 1, 2, 3) being three mass eigenvalues, we can obtain the CKM and MNSP matrices as follows: 

V = {pmHpM 



OiPyO^ , 



U = {PiOi)\P,0, 



ojPuO, , 



(4) 



In this treatment we have assumed massive neutrinos to be the Dirac particles because the overaU neutrino mass 
matrix Af^ is not symmetric. Of course, one may first apply the same treatment to Mq and then invoke the seesaw 
mechanism to produce a Majorana mass matrix for three light neutrinos. 
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where Py = P^P^ = Diag{e*''^i, 6**^2 , 1} and Pjj = PjPy = Diag{e*'^i , 6**^2 , 1} are two diagonal 



phase matrices in a chosen phase convention. 

FoUowing the same model-building strategy as specified in Ref. [12], here we focus on the 



consequences of M^^'^'^ on flavor mixing. We consider the real mass matrix m2^^'^'^ 
where 



/ 









k\ > 










( ° 


—a 












V 








(5) 



with e" and being dimensionless real parameters describing small and asymmetric corrections 
to Ma ■ Treating and as perturbation parameters will technically allow us to perform an 



analytical diagonalization of M( 



(NNI) 



For simplicity, we omit the flavor index a in the subsequent 



discussions. It is easy to exactly diagonalize M^'^'^^) = m(^) in the limit of = ej, = |13] : 



A1A2A3 



(A1-A2 + A3) ' 



Al - ^2 + A3 

- A2 + A3 ; 



(6) 



and 



<n(0) 



(Ai + A2)(Ai-A2 + A3)(A3-Ai) 



O 



(0) 

12 



AiAslAi + Ag) 



(Ai + A2)(Ai -A2 + A3)(A2 + A3 



O 



(0) 
13 



AiA2(A2 - A^) 



Ai)(Ai -A2 + A3)(A2 + A3) 



(0) 

21 



Ai(A3 - A2) 
(Ai + A2)(A3-Ai) 



O. 



(0) 

22 



A2(Ai + A3 



(Ai + A2)(A2 + A3) 



(0) 
23 



Ai)A3 



(A3-Ai)(A2 + A3) 



O 



(0) 
31 



o. 



(0) 
32 



' Ai(A2- 


Ai)(Ai + A3) 


(Ai + A2)(A3- 


-Ai)(Ai-A2 + A3) 


' A2(A2 - 


Ai)(A3-A2) 



o. 



(0) 

33 



V (Ai + A2)(A2 + A3)(Ai-A2 + A3) 

' A3(Ai + A3)(A3-A2) ~ 
(A3-Ai)(A2 + A3)(Ai-A2 + A3) ' 



(7) 
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where the superscript "(0)" imphes the Fritzsch (or e„ = = 0) limit. Note that the above results 
hold for the normal mass hierarchy (i.e., A^^ < < A3) 

Switching on the corrections of M^^^ to M^^\ one may calculate O and O' appearing in the 
bi-unitary transformation M^^^^^O' = M by following a perturbative way. We have done 
such a perturbative calculation both to the first order of e„ and and to the second order of 
them, in order to examine whether the first-order analytical approximations are good enough. Of 
course, the fermion mass hierarchies should also be taken into account in our calculation. Given 



<C TTT-c <C TTij and <C rrig <C m^, it is easy to simplify Eq. (7) by 



making reliable analytical approximations. As only the normal hierarchy of three neutrino masses 
(i.e., nil < 771-2 < "^3) is allowed in this scenario, it is also straightforward to simplify Eq. (7) 
in the neutrino sector. This treatment might not be excellent if three neutrinos have a relatively 
weak mass hierarchy, but it should be good enough for us to reveal the salient features of non- 
symmetric corrections to the Fritzsch textures. Our second-order analytical approximations in 
diagonalizing M^^^^^ = M^^^ + M^''\ which include the 0{e'^) and 0{ef) corrections, support the 
above arguments but they are too complicated to be presented here. To the first order of e^j and 
ej,, we simply take O = O^^^ (1 + X) and O' = O*-*'-* (1 — X) with X being anti-Hermitian (i.e., 
X^ = —X) and proportional to the perturbation parameters and e^,. In this case, 

m(f) = 0(0) M O(0)^ , 

= 0(0) (^XM + MX^ 0(°)^ . (8) 

Then we can determine the matrix elements of X in terms of those of O^^^ and the perturbation 
parameters e„ and e^. After an algebraic calculation, we obtain 



O 



11 



o 



(0) 

11 



o 



12 



o 



(0) 

12 



o 



13 



o 



(0) 
13 



o. 



21 



o 



(0) 

21 



o. 



22 



o 



(0) 

22 



1 + 



1 + 



1 + 



Ai [(Ai-A2)^ + (2Ai-A2)A3 + A^ 
(Ai-A2)(Ai + A3)(Ai-A2 + A3) 

A2 [Af + (A2-A3)^-Ai( 2A2-A3)] 

A2 + A3) 

- As)^] A3 



(Ai-A2)(A2-A3)(Ai- 
[A2-Ai(A2-2A3) + (A. 



(A2-A3)(Ai + A3)(Ai -A2 + A3) 
X2X^{2X-^ - A2 + A3) 



Ai 


Ai 


-A2 


+ A3 




A2 




Ai 


-A2 


+ A3 




A3 





An + Aq 



O. 



23 



o 



(0) 
23 



o. 



31 



o 



o. 



32 



o 



(0) 
31 

(0) 
32 



1 + 



o. 



33 



o 



(0) 

33 



(Ai 


-A2)(Ai + A3)(Ai-A2 
AiA3(A^ - 2A2 + A3) 


+ A3) 


(Ai 


- A2)(A2 - A3)(A]^ - A2 
AiA2(Ai — A2 + 2A3) 


+ A3) 


(Ai 


+ A3)(A2-A3)(Ai-A2 
A2(A2 — A3)A3 


+ A3) 


(Ai 


-A2)(Ai + A3)(Ai-A2 
AiA3(Ai + A3) 


+ A3) 


(Ai 


- A2)(A2 - X^){Xi - A2 
Ai(Ai - A2)A2 


+ A3) 


(A2 


-A3)(Ai + A3)(Ai-A2 


+ A3) 



e« + 



Ai 


Ai 




A2 + A3 






A2 


Ai 




A2 + A3 






A3 


Ai 




A2 + A3 




A2 


-A3 


Ai 




A2 + A3 




Ai 


+ A3 


Ai 




A2 + A3 




Ai 


-A2 


Ai 




A2 + A3 



(9) 



This mass hierarchy is consistent with both the observed mass spectra of charged fermions and the expected 
mass spectrum of neutrinos. Although an inverse mass hierarchy is also possible for three neutrinos, it cannot be 
consistent with the Fritzsch texture M^^' \W\ or its non-Hermitian extension under discussion. 
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Given A^^ ^ A3 and A2 ^ A3, the above expressions of Oj • (for i,j = 1, 2, 3) may approximate to 



Oi3 ^ 



'22 



(1 ~ 



^^^"-\/a;ta^('-a^^"-^^^'' 



^^^^ ^ 1 + T^^a + 



'3 



(Ai + A2)(A3-Ai + A2) V A2-A1 ' A 



A3 ~ A^ + A2 '''' 





Ai(A2- 


-Ai) 


(Ai 


+ A2)(A3 


- Ai - A2) 




A2(A2 - 


-Ai) 


(Ai 


+ A2)(A3 


+ Ai + A2) 



It is obvious that the off-diagonal matrix elements of O are more sensitive to the corrections 
induced by the perturbation parameters e„ and e^. 



3 A numerical illustration 

Now let us take a look at how sensitive the flavor mixing parameters of leptons and quarks are to 
the perturbation parameters and e'^ (for a = u, d; /, v). We first discuss the CKM matrix V and 
then analyze the MNSP matrix [/ in a numerical way. 

(A) The CKM matrix V 

Given six quark masses as the input parameters, the CKM matrix V = OtPyO^ still contains 
six free parameters: e^, e^, e'^, ef, (f>i and ^2- Because the four perturbation parameters must 
be small, we require je^^'Jj'l ^ 0.3 as the reasonable bounds in our numerical calculation. Then 
the experimental data on four independent observable quantities of V ^ typically chosen as \V^g\, 
l^cfel) \^ub\ ^"^^ sin 2/3 with /3 = arg [—(V^^*^^)/{y^^V^^y\ being an inner angle of the CKM unitarity 
triangle [2], will allow us to constrain the parameter space of quark mass matrices Mu'^^^^ and 
^^NNi)^ Such a constraint will be useful for model building. 

To simplify the numerical calculation, we fix the values of quark masses at the electroweak 
scale ^ = Mz as follows: m„ = 2.0 MeV, = 0.557 GeV, = 168.3 GeV; = 2.7 MeV, 
m, = 47 MeV and m^, = 2.92 GeV [H [12]. In addition, we adopt |KJ = 0.2255 ± 0.0019, 
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\VJ = (41.2±1.1) X 10-3, \V^^\ = (3.93 ± 0.36) X 10-3 and sin 2/3 = 0.681 ± 0.025 [2J. By inputting 
the chosen values of six quark masses and allowing six free parameters of V to vary, one may 
then obtain the outputs of {V^J, \V^i,\, and sin 2/3 which are required to lie in their respective 
ranges given above. This treatment leads us to the parameter space of e" versus versus 
and (j)^ versus (j)2, as shown in Fig. 1. Some comments and discussions are in order Q 

(1) Although we have set the bounds |e"'^| < 0.3, their allowed ranges are actually much smaller. 
In particular, e" / and e'^ ^ hold, but ej^ and are possible to vanish. This observation implies 
that both M,!^^^^ and M^^^'^^ must be non-Hermitian. On the other hand, one can see that < 
and > hold. In comparison, is negative in most cases and can be either positive or 
negative. To reduce the number of free parameters from four to two, one may either switch off eJJ 
and or set e" = and = e^, or take eJJ = — and = — e^, and so on. Such assumptions will 
strictly constrain the textures of quark mass matrices and might be suggestive for model building. 

(2) It is impressive that two CP-violating phases (f)i and (j)2 are restricted to a quite narrow 
parameter space. In particular, ~ 90° and \4>2\ ~ 0° imply that (j)^ dominates the strength of 
CP violation in the CKM matrix V. This feature is similar to the one showing up in some Hermitian 
modifications of the Fritzsch ansatz, such as the four-zero textures of quark mass matrices |15j . 

(3) In the quark sector we follow Ref. [12] to define the small parameter 

'^^yie-af + ielf + ieif + ietf (11) 

to measure the overall non-Hermitian departure of Mu^'^^^ and M^^'^^^ from the Fritzsch texture. 
The Jarlskog invariant of CP violation [16], defined as J for the CKM matrix V, can be calculated 
via J = ^^iV^gV^i^V*^V*g) . We illustrate the numerical dependence of J' on e in Fig. 2, where 
J ~ 3.7 X 10-^ for e ~ 0.2. The sign of J is fixed by that of sin 2/3. 

(B) The MNSP matrix U 

Given three charged-lepton masses and two neutrino mass-squared differences Am|]^ = m| — ml 
and Am^2 = "^1 ~ ""^i' 

the MNSP matrix U = OjP^O^ depends on seven free parameters: e^, 
e\j, e'^, e^, m^, tp^ and ip2- Again we require le^j'^^l < 0.3 as the reasonable bounds. The present 
neutrino oscillation data on three flavor mixing angles, denoted as 6i2, ^13 and ^23 iii the standard 
parametrization of U (i.e., tan0]^2 = \Ue2/^ei\^ sin^^^g = \U^^\ and tan023 = I^^3/^t3I) i^J, will 
allow us to constrain the parameter space of lepton mass matrices Adj^^^^^ and mI)^^^\ For 
simplicity, we fix the values of three charged-lepton masses at the electroweak scale as follows: 
mg = 0.48657 MeV, = 102.718 MeV and = 1746.24 MeV [14]. Moreover, we assume 
= 0.0025 eV and take Amj-^ = 8.0 x 10^^ eV^ and Am|2 = 2.5 x 10"^ eV^ together with 
30° < 6-^^2 < 38°, 36° < 6*23 < 54° and < 10° in our numerical calculation. Then the normal but 
weak neutrino mass hierarchy is measured by two mass ratios mi/m2 — 0.27 and ra2/ra^ ~ 0.18. 
By inputting the chosen values of charged-lepton and neutrinos masses and allowing the unknown 
parameters of U to vary, one may obtain the outputs of 0^2; ^13 ^'^^ ^23 which are required to lie 
in their respective ranges given above. This treatment leads us to the parameter space of e\ versus 
e^, versus and versus ^P2i shown in Fig. 3. Some comments and discussions are in order. 



! have recently analyzed the CKM matrix V in such a way [12] . Our more comprehensive 
analysis not only confirms their results but also provides ourselves with a meaningful calibration as we extend the 
same analysis to the lepton sector. Our results for the MNSP matrix U in section 3 (B) are completely new. 
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(1) Because of m^/m^ <C m-^^/m2 and m^/m^ <^ 1712/171^, the MNSP matrix U is expected to 
receive more contributions from the neutrino sector rather than the charged-lepton sector. That 
is why the bounds on |e^^| are much looser than those on |ej^^|, as one can see from Fig. 3. We 
find that e'^ and e'j^ are negative in most cases. To reduce the number of free parameters from four 
to two, one may simply switch off and e\^. Although it is also possible to set = £5 = 0, it 
will be impossible to get ^23 — ^^5° for the atmospheric neutrino mixing angle [10]. More precise 
data to be extracted from the upcoming neutrino oscillation experiments may help us constrain 
the ranges of ^ and ^ to a better degree of accuracy. 

(2) Current neutrino oscillation experiments set no constraint on the CP-violating phase ipi- 
But the other CP-violating phase (^2 is well restricted to be around 180°, as shown in Fig. 3. The 
reason is simply that a sufficiently large value of ^23 requires ~ 180°. To see this point more 
clearly, we write out 



tan 623 







Ur3 





(12) 



;i + h^ei + hlel) + h,K (1 - 4) (1 - e''^^ 

with the help of Eq. (10), where /i^ = {m^ — m^)/m^ ~ ^m^/m^ and /i^ = (rrig — m^)/m.^ . 
It becomes transparent that ip2 ~ 180°, together with e^, < and < 0, may enhance the 
magnitude of tan ^23 ^^id make ^33 closer to its best-fit value (i.e., 633 — 45°). 
(3) Similar to the definition of e in the quark sector, a small parameter 



^' - \ ^{^if + {Af + {^'if + {^lf (13) 

can also be defined to measure the overall non-Hermitian departure of M^^^'^^ and M^'^^^^ from the 
Fritzsch texture. We illustrate the dependence of three flavor mixing angles on e' in Fig. 4. The 
Jarlskog invariant of leptonic CP violation, which can be calculated via J' = Im(C/g2^/x3 ^1=3^/^2)' 
also shown in Fig. 4. We see that it is possible to reach \J''\ 1.5 x 10^^. The CP-violating effects 
at this level should be observable in the future long-baseline neutrino oscillation experiments. 

Finally, it makes sense to compare between the MNSP and CKM matrices derived from the 
same NNI textures of lepton and quark mass matrices. Given |e"'^| < 0.3 and \e^'^fj\ < 0.3, it 
can be concluded that the smallness of three quark mixing angles is primarily attributed to the 
strong mass hierarchies of up- and down-type quarks, while the largeness of solar and atmospheric 
neutrino mixing angles is mainly ascribed to the relatively weak hierarchy of three neutrino masses. 
Of course, the CP-violating phases play an important role in either the lepton sector or the quark 
sector. Note that 1/^^ is smaller in magnitude than all the other elements of the CKM matrix V, 
and C/g3 is also the smallest element of the MNSP matrix U. In other words, ^^^3 is the smallest 
mixing angle in both lepton and quark sectors. This interesting feature is a natural consequence 
of the flavor textures of leptons and quarks together with their corresponding mass hierarchies. In 
addition, the fact that |e"'^| ~ je^^^^l ~ 0.2 is favored by current experimental data should be quite 
suggestive of a unified flavor model of leptons and quarks. 



4 Summary 

We have introduced non-Hermitian and nearest-neighbor-interacting perturbations to the Fritzsch 
textures of lepton and quark mass matrices such that both of them can fit current experimental 
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data very well. In particular, we find that it is possible to obtain — 45° for the atmospheric 
neutrino mixing angle and predict 9^^ ~ 3° to 6° for the smallest neutrino mixing angle when the 
dimensionless perturbation parameters in the lepton sector are at the 20% level. We have shown 
that the same level of perturbations is required in the quark sector, where the Jarlskog invariant 
of CP violation is about 3.7 x 10~^. In comparison, the strength of leptonic CP violation is likely 
to reach about 1.5 x 10"^ in neutrino oscillations. 

As shown in Ref. [12] . the NNI texture of quark mass matrices can be derived from the 
introduction of an Abelian flavor symmetry (e.g., the minimal realization of this idea requires a 
flavor symmetry in the context of a two-Higgs doublet model). We can follow the same procedure 
to obtain the NNI texture of lepton mass matrices if massive neutrinos are the Dirac particles. In 
the presence of a few heavy Majorana neutrinos, one may first impose the aforementioned flavor 
symmetry on the Yukawa interaction of neutrinos to get the NNI texture for the Dirac neutrino 
mass matrix and then achieve the Majorana mass matrix for three light neutrinos via the seesaw 
mechanism. There are therefore a number of possibilities of model building, but the numerical 
results must be very different from what we have presented in this work. One may explore such 
possibilities once more accurate experimental data on neutrino masses and lepton flavor mixing 
angles are available in the (near) future, and in particular when the simple scenario discussed in 
this paper is phenomenologically discarded or becomes less favored. 

Let us reiterate that a parallel study of lepton and quark mass matrices, such as the approximate 
Fritzsch textures under discussion, is useful to reveal the similarities and differences between the 
lepton and quark sectors. It should also be helpful for building a unified flavor model of leptons 
and quarks with the help of proper flavor symmetries. 

This work was supported in part by the National Natural Science Foundation of China under 
grant No. 10875131 and in part by the Ministry of Science and Technology of China under grant 
No. 2009CB825207. 
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Figure 1: An illustration of the parameter space of e" versus ej^, versus and (p-^ versus 
constrained by current data in the quark sector. 
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Figure 2: An illustration of the dependence of the Jarlskog invariant of CP violation J on 
overall perturbation parameter e in the quark sector. 
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Figure 4: An illustration of the dependence of three flavor mixing angles and the Jarlskog invariant 
of CP violation J'' on the overall perturbation parameter e' in the lepton sector. 
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